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nom de plume Student. The approach that we present can be used to price European options 
using other distributions and yields the Black-Scholes formula for returns described by a normal pdf. 
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I. INTRODUCTION 



<N 

In the fall of 2008, many investors and analysts witnessed multiple "once-in- a- lifetime" events in a single week, with 
disastrous consequences in many cases for the financial health of their portfolios. The once-in-a-lifetime designation 
c*j . of the events arises from a calculation of the probabilities of events based on normal statistics. The returns were more 
than several sigma beyond the expected return, and by normal statistics, thought to be impossible. 

The difficulty lies in the use of a normal pdf to describe returns. It is known that returns have "fat tails" [J0, but 
for mathematical convenience and perhaps force of habit, a normal pdf is often applied. The celebrated Black-Scholes 
formula for pricing European options is based on several assumptions, one of which is that the returns are described 
by Brownian motion 0, |j, Q . The underlying pdf for Brownian motion is a normal pdf. 

In this paper we price European options using a log Student's i-distribution. We call these formulae of the prices 
of European options Gosset formulae, in honour of the author behind the nom de plume Student. We present 
evidence that stock returns are fit by a Student's ^-distribution, in agreement with known results [1, 0, IS 0] ■ We 
then demonstrate the "fat tails" of the i-distribution. One of the difficulties with the fat tails is that one of the 
integrals, which is required to price an option, diverges. We present two similar approaches to handle successfully the 
divergence, and hence find a price for the option. These approaches can be used to price options using pdf's other 
than the i-distribution. However, we restrict our attention to the ^-distribution and compare prices for options using 
the Gosset and Black-Scholes formulae. In general, the Gosset formula yields prices for options that are higher than 
• • ■ the prices found from the Black-Scholes formula. This is not unexpected, as the Gosset formula gives greater weight 
. £h t° events in the tails. 

The starting point for the pricing of the European options is the arbitrage theorem [l(J EH E3] • 



II. ARBITRAGE THEOREM 



Consider an experiment that returns riie^) for a bet of unity on outcome ei of the N possible outcomes ej of the 
experiment. In short, the arbitrage theorem states that either the mean return as calculated over the probability of 
the events ej, E{n(ej)}, is zero or there is a betting scheme that leads to a sure win. The betting scheme may allow 
for negative, zero, and positive bets. Puts and calls allow for negative, zero, and positive bets. 

Let ej be N disjoint events of which one and only one must occur. Let pi be the probability of occurrence of event 
et and c = [c\, Ca, ■■■Cn) be a vector of bets on the N outcomes. The bets Cj may be zero, positive, or negative. The 
gain for a wager of Cj on the i th outcome given that event ej occurred is gi(ej) — Cj x rj(ej). Note that in a typical 
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bet, ri(ej) = — 1 if i ^ j; the player forfeits the wager a if event i does not occur. The mean gain for a wager on the 
jth even t i s g i — Ci x fi and is the expectation over the ej, which is obtained as ^2jPj X 9i( e i) = c i x ^2jPj x r i( e j)- 

Following the example of de Finetti [l(| , we recognize the <ji as N equations in N unknowns Cj . To ensure that no 
arbitrage exists, the average gain g~i = 0. This is only possible for arbitrary Cj if fj = 0. If f,; ^ 0, then the Cj can be 
selected to yield, on average, cji > 0. Thus to ensure no sure win, fj = B{rj(ej)} = 0. That this condition must hold 
for all i is a statement that there are no lucky numbers or values for assets. In the context of pricing an option on an 
asset, the events e, could be the final price of the asset falling in a range. 

Ross 12j shows how to price an option using the arbitrage theorem, assuming that the theorem can be generalized 
to handle functions rather than N discrete outcomes. If the cost of an option is C, then by the principle of no-arbitrage 
and the arbitrage theorem C — E{pi}, where the return for a European call option is pi — (Si — K, 0) + and the return 
for a European put option is pi = (K — Si, 0) + . The notation (K — S, 0) + means the maximum value of (K — S) and 
0, with K being the strike price and S the value of the asset. 

Ross [HJ obtains the Black-Scholes formula by assuming log-normal statistics for the value of the asset and making 
some assumptions regarding the time value of the option. It is not necessary to solve the Black-Scholes partial 
differential equation to obtain the Black-Scholes formula for the price of an European option [l2|, EH ■ 

We accept as a starting point that the arbitrage theorem can be generalized from N discrete outcomes to hold for 
functions. This we consider a reasonable starting point. The probability functions and values of assets can always be 
quantized to have N discrete ranges. We price European options assuming log Student's ^-distributions, the arbitrage 
theorem, and results from martingale theory. The results from the martingale theory allow determination of the risk 
neutral price of the option. 

First we demonstrate that the log Student's t-distribution is a reasonable distribution to use to price an option. 



III. STUDENT'S t-DISTRIBUTION 



The "fat-tail" property was one of the earliest observations of asset returns, with published reports dating back to 
the early 1960's [J When the Black-Scholes formula was developed in the early 1970's the shortcomings of the 
underlying normal assumption were identified and attempts were made to incorporate kurtosis, see for example [14j | 
and references therein. Also in the early 1970's the Student's i-distribution was first used to describe asset returns [g]. 
There has been considerable work since 1974 comparing the i-distribution and other distributions to empirical data to 
determine which distribution best fits observations, see for example @, and references therein. The i-distribution has 
been found to provide realistic fits to observed asset returns (particularly equities and equity indices) and in many 
cases provides better fits than other distributions [9(. 

In Fig. 1 we show fits to the observed frequency of occurence of daily return data for the Dow Jones Industrial 
Average (DJIA) (left) and the S&P 500 (right) equity indices. Each plot shows the data with a best fit i-distribution 
and normal distribution. The £-distribution fits the observed distribution well, particularly in the tails. We have 
observed similar quality fits for the i-distribution to frequency of occurence of daily returns for individual equities. 

The fit paramaters with the associated uncertainties, some critical values, and sample descriptive statistics, such 
as the count of data points N, sample mean f, and sample standard deviation s are given in Table I. The critical 
values x c (q) are the returns which solve P{£ < x c } = q. Given that there are 21,186 data points for the DJIA, it is 
expected, using the best fit t-distribution, that only 2 returns will fall below -0.26. This is in agreement with the data 
shown in Fig. 1. 



TABLE I: Fit parameters, critical parameters, and simple descriptive statistics 







DJIA 




S&P 500 




parameter 


normal 




t 


normal 


t 



location p (1.679±1.6) x 10~ 4 (4.023±1.1) x 10" 4 (2.660±1.5) x 10~ 4 (4.270±1.2) x 10~ 4 

scale a (1.162±0.011) x 10 -2 (1.157±0.011) x 10~ 2 (9.55±0.11) x 10~ 3 (8.71±0.19) x 10~ 3 

shape v 2.65 ±0.11 3.24 ±0.19 

|a; c (10- 4 )| 0.043 0.26 0.036 0.19 

|x c (10" 3 )| 0.036 0.12 0.030 0.089 

|a; c (10- 2 )| 0.027 0.053 0.022 0.040 

count 20, 186 14, 870 

f 1.679 x 10" 4 2.660 x 10 -4 

s 0.01157 0.00951 



Figure 2 is a plot of the value of the Student's ^-distribution and the normal distribution as a function of the number 
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FIG. 1: Fit of Student's t-distributions and normal distributions to returns for the DJIA and S&P 500. 

of standard deviations from the true mean. The fat tails of the Student's £-distribution are evident. Events that are 10 
standard deviations from the true value are expected approximately 10 _2 % of the time for the Student's ^-distribution 
with a small number of degrees of freedom, but only 10 _18 % of the time for a normal pdf. For comparison, the values 
for iV(0, 1.5a), a normal pdf with a mean of zero and a standard deviation 1.5x greater than for the other curves is 
presented. A volatility of 1.5x the measured volatility is sometimes used to price options. The enhanced volatility 
effectively gives the normal pdf "fatter tails". It can be observed in Fig. 2 that the N(Q,1.5a) curve tracks the 
Student's £-distribution over a greater distance from the true mean. 

One possible justification for the use of Student's t-distribution to describe the daily returns for an asset is as 
follows. If the returns are assumed to be normally distributed with a given short term variance, the variance is taken 
to be a random variable that varies slowly in time, and ft(x) is the pdf for the return at time some time t in the 
future, then ft(x) — J f t (x\a t ) x pdf at (a t )dat where f t (x\a) is the pdf for the return given a and pdf at (a t ) is the pdf 
for a at time t. Fits to the 22-day volatilities for the DJIA and for the S&P 500 show that a gamma distribution fits 
the 22 day volatilities. A gamma distribution with shape parameter v is chi-squared with 1v degrees of freedom [po] ]. 
and a normal distribution with a variance that is chi-squared distributed is a Student's i-distribution. A Student's 
i-distribution for the pdf of the daily returns is consistent with the data shown in Fig. 1. 

The variance of a Student's i-distribution is vj(v — 2) [15j. As v tends to infinity, the ^-distribution tends to a 
standard normal pdf, which has a variance of one. The v j(y — 2) enhancement accounts for the range of values that 
the variance might take. If it is known that the variance does not change in time (i.e, pdf at (at) = 8(a t ~ &o) where 
8(x) is a Dirac delta function), then the Gosset formula presented in this paper reduces to the Black-Scholes formula. 
In this respect the Gosset formula extends the Black-Scholes formula by removing the condition that the variance is 
constant in time. Under this interpretation, the shape parameter v is a measure of the belief that the volatility will 
not change. A large v means that one believes that the volatility will not change in time. A small v means that 
one believes that there is a broad range of values that the volatility might take. If the volatility is determined as 
the standard deviation of N independent returns that are normally distributed with an unknown variance, then the 
maximum value of v should equal N — 1, which is the number of degrees of freedom in a determination of a sample 
standard deviation. This maximum value of v = N — 1 takes into account the fact that the sample standard deviation 
is a random variable and has an uncertainty associated with it. 

Figure 3 is a plot of the relative areas in the tails for the same distributions presented in Fig. 2. The curves are 
normalised by the area in the tail of a normal pdf. Figure 3 gives the probability of an event x for x > t, relative to 
the probability of the event x for x > t for a normal pdf. The probability of an event lying > 10 standard deviations 
from the mean is 10 20 time mores likely for a process that follows a Student's i-distribution with 3 and 5 degrees of 
freedom than for a process that follows a normal pdf. 
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FIG. 2: Comparison of Student's ^distributions and nor- 
mal pdfs. 
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FIG. 3: Normalised areas in the tails for Student's t- 
distributions with 3 and 5 degrees of freedom and for a 
normal distribution with 1.5x the volatility, each nor- 
malised by the area in the tail of a normal pdf. 



IV. PRICING THE OPTION 



Le St be the price of a stock at time t, t > 0. Let be the strike price at time T, where T is the time when the 
option expires. 

Let St — A t exp(cr t £) be the value of the stock where £ is a random variable. 

The cost of a European call option, calculated at the time of expiration T, is Cr = E{(St ~ Kt) + } = J pdf (St) x 
(St — Kt) + dST , which is the expectation of the maximum value of {St — Kt, 0}. The desired quantity is Co, which 
is the value of the option at time t = 0. The desired quantity can be obtained from the expected time value of money. 
If r(t) is the risk free rate, then Co = E{Ct x exp(J r(t) x tdt)} — Ct x exp(— r x T) when the risk free rate is 
assumed to be time independent. This is a standard assumption in the derivation of the Black-Scholes formula 0,0). 

If £ is normally distributed, then St follows a log normal distribution and the price for the option follows the 
Black-Scholes formula [4j, |fj[. 

If £ follows a Student's t-distribution, then the pdf for £ for v degrees of freedom is 

V df(i)di = i££L x d -L_ = A(„) x d ±-^ , (1) 



r(|) 

which can be written in terms of St — A t exp(cr t £) as 



(' + ?)■ 



A(f) dSt ,„> 

pdf(S t )dS t = — — x *— . (2) 

— at x bt 



(\n{StJA t ))< 



The following integrals are needed to evaluate the price of a call option at the time of expiration, Ct '■ 



/•OO poo 

/ pdf (S T ) x K T xdS T = \ 



K(v) x K 7 



(\n(S T /A T )y 



dST 
<jt x St 



(3) 



in(K T /A T ) 



A ^ X ^ x dZ 



■5 




oc 



A(z/) x K T 



dS T 



/ pdf(S T ) xS T x dS T 



( 1+ (\n(S T /A T 




X 



<7t x 



(4) 



OG 



A T exp(<7 T £) 



\n{K T /A T ) 



AM 



(l + ?r 



X d£ • 



Note that At is a function of So, r, T, x c (p), and Z, and that these variables are not functions of the random variable 
£ (see Eq. (6) and the accompanying text). 

The integral in Eq.(4), which is for a call option, will be infinite for small v\ the exponential term in the numerator 
will dominate the denominator. Viable and physical results are obtained by restricting the upper limit of integration 
(i.e., truncate the pdf) or by restricting the value of the asset St (i.e., cap the value). It is not physical for the value 
of the asset St to approach infinity. 

For a put option this problem does not directly arise as < St < Kt and thus the integrals are evaluated from 
— oo to \xy{Kt / At) / (Jt ■ However, the average price of the asset is needed to price an option, and thus the pricing of 
a put option suffers also from an integral that diverges. 

The integral that diverges is Eq. (4). If the value of the asset is restricted to be < At exp(<jTX c ), then the offending 
integral, Eq. (4), can be rewritten as 



The integral in the last line of Eq.(5) is the probability q = (1 — p) that the stock price St — At exp(aT^) is greater 
than some critical value x c (p), which we abbreviate to x c . The last line in Eq. (5) gives the expected value of the 
cost of capping the value of the asset to At cxp(<ttx c ). 

In pricing a call option, one could also truncate the pdf by setting the last line of Eq. (5) to zero. That is, one 
assumes that the pdf is zero for values greater than the critical value x c . Truncation sets an upper limit on the 
value of the stock: St < x c . It is necessary to divide the truncated pdf by p to maintain the normalisation of the 
truncated pdf. This follows from the definition of conditional probability: P{A n B} = P{A n B}/P{B}. In this 
case, B := < St < S maXl A := St, An B = St for < St < S ma x and equals the null set otherwise. Thus 
P{S T \0 <S t < S max } = P{S T }/P{0 <S t < S max } = P{S T }/p if < S T < S max and = otherwise. 

Table II gives critical values and the maximum increase in stock price = cxp(cttx c ) for a normal pdf and a Student's 
^-distribution with v = 5 degrees of freedom, volatility &t — 0.40, and various value of p. 

The effect of the fat tails of the Student's ^distribution are evident. For a confidence level of 99.99%, the critical 
value for a Student's f-distribution with 5 degrees of freedom is 9.678 and the stock price has increased at most 
47.992x, a 4 798% increase. For a normal pdf, the critical value is 3.719 and there is a 99.99% chance that the stock 
price has increased at most 4.426x, a 442% increase. The other interpretation is that there is only a 0.01% chance 
that the stock price has increased by > 47.992x under the Student's ^-distribution or > 4.426x under the normal 
distribution. 

A confidence level of p = 0.9999 (i.e., 99.99%) seems unrealistic. This confidence level is one part in 10,000 and 
allows for increases of cxp(<ttx x ) in the price of the asset that may be physically impossible to achieve (c.f. Tabic 
III). If the time period is one year, then once in 10,000 years is too large an extrapolation to be meaningful. Even 
p = 0.999 seems unrealistic. The best way to price options may be to pick an upper bound x c = \n(S max )/<TT (or 
confidence level p) for the price of the asset based on insight and agreement between parties, and then calculate the 
price using a capped distribution. 

Although a confidence level p — 0.9999 might be unrealistic or unphysical in terms of expected increase in the 
stock price and extrapolation, the results presented here show that the methods allow a price for the option to be 
calculated. The maximum increase in value of the asset is a function of the number of degrees of freedom and the 
assumed confidence level p or critical value x c (p) . Table III presents the critical values for different values of p and 
degrees of freedom v. Note that a Student's ^-distribution with v = oo yields a normal pdf and that a ^distribution 
with v > 40 is often considered to be a normal pdf. 
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The variable At is required to price the option. This value is obtained from the assumption that an option is a fair 
wager (i.e., a martingale), as shown next. 



TABLE II: S t /A t at x c (p) for volatility = 0.4 and selected p. 





Student's t, v = 5 


normal 


p 


x c (p) 


exp(a T x c ) 


%c(p) 


exp(cr T x c ) 


0.9 


1.476 


1.805 


1.282 


1.670 


0.95 


2.015 


2.239 


1.645 


1.931 


0.99 


3.365 


3.842 


2.326 


2.536 


0.995 


4.032 


5.018 


2.576 


2.801 


0.999 


5.893 


10.56 


3.090 


3.442 


0.9999 


9.678 


47.99 


3.719 


4.428 



TABLE III: Critical values x c (p) for various values of p and v. 



V 


p = 0.90 


p = 0.95 


p = 0.99 


p = 0.999 


p = 0.9999 


3 


1.638 


2.353 


4.541 


10.21 


22.20 


4 


1.533 


2.132 


3.747 


7.173 


13.03 


6 


1.440 


1.943 


3.143 


5.208 


8.025 


40 


1.303 


1.684 


2.423 


3.307 


4.094 


CO 


1.282 


1.645 


2.326 


3.090 


3.719 



V. MARTINGALES, DOOB DECOMPOSITION, AND RISK NEUTRAL PRICING 

A martingale is a mathematical model for a fair bet[IB|. A fair bet can be defined as one for which the conditional 
expectation (i.e., given all the knowledge at the time of the bet) of the payout of the bet should be zero. 

A martingale is defined such that E{X n+ i\F n } = X n for all n where F n is the knowledge (history) up to time n and 
the knowledge increases as n increases. The X n form a stochastic process. A martingale is non-anticipating because 
the outcome of a martingale depends only on the past and not on the future. 

A submartingale is defined such that X n < E{X n+ i\F n } and a supermartingale is defined such that X n > 
E{X n+ i\F n }. A submartingale represents a favourable process (the expected outcome is favourable) and a super- 
matingale represents an unfavourable process (like a betting at a casino). Casinos offer submartingales (the patrons 
play supermartingales) and financiers aspire to offer or to play submartingales. 

If X n is a martingale, then exp(X n ), exp(— X n ), \X n \, and X% are all submartingales. 

The Doob decomposition shows that any submartingale X n is the sum of a martingale M n plus an increasing 
process, X n = M n + A n , where A n < A n+ i and A n depends on the knowledge to n — 1, F n -i 16]. 

Let St be the spot price at time t for a stock. St = So exp(a x t) where a is the expected return on a non-dividend 
paying stock, a is from a capital asset pricing model (CAPM) or equivalent p|. (a — r) is called the risk premium and 
a = r + (a — r) is the expected return on the stock, which equals the compensation for the time value of the money 
invested in the stock plus the compensation for the risk of holding the stock. The time value is separated out since 
the value of the option Ct is calculated at time T, but the value at time 0, Co, is desired, and Co = Ct X exp(— r x T), 
assuming the risk free rate r is known and constant in time. 

The Doob decomposition is used to create a martingale for pricing options on the stock. The drift of the stock 
is subtracted from the process. This procedure of subtracting the drift from the process also leads to risk neutral 
pricing, where one is able to price an option without explicit knowledge of the risk premium. 

The average value of a stock at time t = when we wish to price an option is Sq. The average price of the stock 
at some time later is E{St} — E{A t exp(<7t£)} where £ is a random variable and the drift is contained in A t . For 
a martingale (a fair wager), E{St} = So x exp(r x t). One may use the Doob decomposition to offset the drift and 
require E{A t exp(crt£)} = {A t exp(cr t (£ — £o))} = So x exp(r x t). This is equivalent to requiring that the pdf for St 
should be centred about the drift owing to the risk premium. The exp(r x t) term takes into account the time value 
of money and allows for comparison of values at two different points in time. 

It is perhaps simpler to find the normalization A t exp(— <7t£o) that satisfies the martingale condition than to shift 
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the pdf. For a capped distribution, 



pdf(S T ) xS T xdS T = I A(u) x A t^pM^-Jo)) x ^ (6) 



> A T x / AM x ""'^it x ,,( 



do 



1 + V 



At x exp(<7T# c ) x / A(^) 



1 + V 



f+1 

t?\ 2 



Let Z be the value of the integral in the middle row of Eq. (6) and let (1 — p) be the value of the integral of the last 
line of Eq. (6), i.e., (1 — p) = P{£ > x c }. Equation (6) can then be rewritten as 

S x exp(r xT) = A T x Z + A T x exp(cr T x c ) x (1 - p) , (7) 

and solved for the value of At'- 

A T = 5oXe " p(rXr) , . (8) 
Z + (1 —p) x exp((7tx c ) 

For a truncated distribution where the pdf is set = for £ > x c , At = So x exp(r x T)/Z and the pdf in the integral 
that defines Z (middle line of Eq. (6)) will have a factor of 1/p where p — P{£ < x c } is the probability that £ < x c . 
This follows from the definition of conditional probability: P{x < £ < x + dx\£ < x c } = P{x < £ < x + dx H £ < 
X c}/P{£, < %c} — pdf(£)/p if £ < x c and zero if £ > x c . 

Note that Z = Z(a, x c (p), t, v) and that Z is not a function of the random variable £. 

For a log normal distribution with fj, = and a = 1, x c = +oo,p = 1, Z = exp(a T /2), and At = So X cxp(r x T — 
a T /2). Thus C T = S x exp(r x T) x N({\n{K T /S ) - r x T- a T /2)/a T ) - K T x N{{\n(K T /So) ~ r x T + ct t /2)/ct t ) 
where i\T(x) is the cumulative distribution function (CDF) for the unit normal pdf, N(x) = P{X < x\a — 1}. 

Co = C T x exp(-r x T) = 5 x iV(d x ) - exp(-r xT)xKx N{d 2 ) where d x = (bx(S /K T ) + rxT + a T /2)/a T 
and d 2 = d\ — ot- This is the standard form for the Black-Scholes formula for a European call [3, Q. Note that 
Ko = ifTexp(— r x T) and the expressions can be rewritten in terms of Ko : Co = So x N(d\) — Kq x N(d2) with 
rfi = (ln(5o/ifo) +o't/2)/ 't and d 2 = di — a T - 

Thus the value for appears correct and it should be possible to price an option assuming a Student's ^-distribution 
(or any other distribution) using the approach and equations developed in this manuscript. 

The price at the time of expiration for a European call, Ct, using a Student's f-distribution and capping the value 
of the asset at a critical value x c is 

C T = A(i/) x / — j x d£ (9) 

del 



v 1 " / 

+ A(u) x (A T exp((TTX c ) — K T ) x 



X < (l + f)' 



For a capped asset, Ay is defined by Eq. (8). 

For a truncated i-distribution, the price of a call Ct is obtained by dropping the last line of Eq. (9), dividing 
the integral in the first line of Eq. (9) by p(x c ) — P{£ < x c }, and using At = So x exp(r x T)/Z. It is helpful to 
remember that Z is defined by some of the integral on the first line of Eq. (9) (see Eqs (7) and (6) for the definition 
of Z) and thus is also divided by a factor of p(x c ) for a truncated pdf. The value of the integral in the last line of Eq. 
(9) equals p(x c ), which we shorten to p. 

The price for a European put option at time T is Pt and equals 

ln(K T /A T ) 

... f "t Kt — At exp((Tr£) , „, 

P T = K{v) x / -1 1 — x dC , (10) 
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for both a capped asset and a truncated pdf, where the truncation occurs at x c > \u(Kt / At) / <tt ■ The value of Z is 
different for the capped and truncated approaches, and this makes At different for the two approaches. 

In pricing a put, one could also truncate the distribution (i.e., set the lower limit of integration to a critical value 
x p and set the pdf to zero for values < x p ) or use a floor at the critical value (rather than a ceiling as was done for a 
call). The values of Z and At would be modified in a straightforward manner to accomodate the floor or truncation. 
In this paper we have chosen to use the full range of the integral in Eq. (10). There are six possibilities to calculate an 
option: two choices for the call (cap or truncate — one of which must be chosen) and three choices for the put (floor, 
truncate, or nothing). Our goal is to present the approach to pricing an option using a Student's i-distribution. In 
this regard, is it necessary to truncate or to cap the distribution at the upper end. We leave discussion of the merits 
of a floor and truncation to future work. 

We consider a formula for the price of an European option that is obtained using either a maximum value for the 
asset or a truncated Student's ^distribution to be a Gosset formula. 

There is little difference in philosophy between options priced with a capped value or with a truncated pdf. In 
both approaches the return is assumed to be always less than a critical value. With the capped asset approach, more 
weight is applied to the maximum value of the asset than for the truncated pdf approach. The truncated pdf is 
more compact to write as there is not a second, additive term. There is, however, a difference in risk between the 
two approaches. The truncated pdf approach assumes that the value of the asset will never exceed the critical value. 
There is no guarantee that the value will not exceed the critical value. 

The approach of capping the payout of a benefit is common in insurance. Capping the payout reduces the cost to 
purchase the benefit. Capping the payout removes the risk associated with the benefit exceeding the critical value. It 
is known, by design, that the payout will never exceed the cap but it is accepted that there is a non-zero probability 
that the expense will exceed the benefit. 

Capping the value of the asset underlying a call option removes the risk associated with the value of the asset 
exceeding a critical value (assuming that the buyer of the call option agrees to the condition that the payout not 
exceed the cap value), and thus reduces the cost of the option. Capping the value is one approach to keep the integrals, 
which are needed to price an option, finite. The capped Gosset formula would be the proper formula to use to price 
a European option for an asset that is capped in value. 

Put-call parity holds regardless of the underlying pdf. Consider E{St — Kt} = E{(St — Kt) + — (Kt — St) + } = 
Ct — Pt where (x) + — max{x, 0}, Ct is the cash value of a call evaluated at time T, and Pt is the cash value of a 
put evaluated at time T. 

For a martingale (a fair wager), E{S t } — S x exp(r x t), and thus E{St - K T } = S x cxp(r x T) — K T = Ct — Pt- 
From the time value of money, x t = x x exp(r x t), then Co — Po = So — Kt x exp(— r x T) or Co — Pq = So — K . 

One of the confusing features of pricing an option are the time values. Typically the option is priced at time = 
for a stock price at time = and a strike price at time = T. One is usually given S, a, and r implicitly for time = 
(i.e., So, and r = r(0)), the strike price K implicitly for time T, (i.e., K = Kt ), and asked implicitly for C and P at 
time = (i.e., Co and P ). The volatility a is evaluated at time = 0, and is propagated to time = T using diffusion: 
<jt — f x VT. For diffusion, the variance increases linearly with time, ot is used in the pricing of the options. For 
pricing an option with the Gosset formula, it is necessary to specify the shape parameter v, which is a statement of 
belief about the range of values that the volatility might take at time = T. 

VI. COMPARISON OF PRICES FROM GOSSET AND BLACK-SCHOLES FORMULAE 

Figure 4 is comprised of plots of the price of a call option calculated using the Gosset formula minus the price 
of the same option calculated using the Black-Scholes formula, as a function of the number of degrees of freedom 
assumed for the Student's t-distribution. A capped distribution was used for the Gosset formula and confidence levels 
p of 0.99, 0.999, and 0.9999 were used. For small numbers of the degrees of freedom v and p — 0.9999, the Gosset 
formula predicts a price difference > $5 between the values predicted by the Gosset and Black-Scholes formulae. For 
the calculations, a risk-free rate r = 3%, Sq = $50.00, Kt = $49.00, ot = 0.3, and T = 1.0 year were assumed. The 
Black-Scholes formula predicts a price of $7.12 for a call option with these parameters. The Student's ^-distribution 
is considered to be equivalent to the normal distribution for the number of degrees of freedom v > 40. The differences 
between the prices obtained from the Gosset formula and the Black-Scholes formula are between $0.06 and $0.11 for 
40 degrees of freedom for confidence levels of 0.99, 0.999, and 0.9999. 

Figure 5 is comprised of similar plots as for Fig. 4, except that a truncated distribution was used to calculate the 
price of the call option with the Gosset formula. The price of a call option obtained by using a truncated distribution 
is less than the price obtained by using a capped distribution. For a confidence level of 99% and v > 25, the 
Gosset formula predicts a value for a call option that is less than the value predicted by the Black-Scholes formula. 
This is not surprising given that for large v, the Student's ^-distribution approaches the normal pdf. As the two 
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distributions become similar, any approach that excludes data in the tails will predict a lower value for the option. 
The difference in cost of the option between the capped asset approach and the truncated distribution approach owes 




FIG. 4: Price difference, Gosset — Black-Sholes, as a func- FIG. 5: Price difference, Gosset — Black-Sholes, as a func- 
tion of v for confidence levels of 0.99, 0.999, and 0.9999 tion of v for confidence levels of 0.99, 0.999, and 0.9999 
using a capped distribution. using a truncated distribution. 



to contributions from three terms. The Z term is different for the two approaches and this leads to differences in the 
integrands and lower limits of integration. For p = 0.9999, v = 3, r = 0.03, or = 0.3, T = 1.0, Sq = $50.00, and 
Kt = $49.00, the price difference, capped minus truncated, is $1.48. The contribution to cost from the values above 
the cap (i.e., the contribution from the last line of Eq. (9)) contributes $3.14 to the difference. The integrand for 
the truncated distribution is larger than for the capped asset. This contributes — $1.52 to the difference. The lower 
limits of integration are slightly different for the capped asset (Z — 1.281 yielding \vl[Kt I At) I <Jt = 0.6583) and for 
the truncated distribution (Z = 1.203 yielding \ii(Kt / At) / <Jt = 0.4488). This area contributes —$0.10, for a total 
difference of $1.52 in the costs of the call option when calculated using the two approaches. 

Figure 6 gives the value of a call option using the capped Gosset formula minus the value of a call option computed 
using the Black-Scholes formula, as a function of the confidence level. For large numbers of degrees of freedom v : the 
Gosset formula predicts a value that is less than the value predicted by the Black-Scholes formula. Again, this is not 
surprising as the Student's i-distribution will approach the normal pdf for large v and any pricing mechanism that 
excludes area in the tails will predict a lower price. 

Figure 6 demonstrates that the price of an option for a ^-distribution with a small v increases dramatically as the 
confidence level increases. For a confidence level approaching 100%, the price tends to infinity since one of the integrals 
required to price the option diverges. Although we use the term confidence level p, it is helpful to remember that the 
confidence level gives a critical value, x c (p). The exponential of the critical value, exp(aTX c ), gives the value that the 
asset will not exceed in a capped asset approach, or is hoped will not be exceeded in a truncated pdf approach. 

Figure 7 presents similar information as Fig. 6 and is presented for completeness. The data for Fig. 7 was calculated 
using a truncated pdf rather than a capped value for the asset. 

Figure 8 gives the prices of calls and puts as a function of the t = price of the asset. The capped Gosset formulae 
for Ct and Pt, Eqs. (9) and (10), were used to calculate the prices of the options. When the option prices were 
used in the put-call parity formula, the difference was of order of the precision of the calculation. This held for both 
the truncated and capped approaches. Figure 8 displays prices for v = 3, v = 5, and for the Black-Scholes formula. 
A confidence level of 99.9%, K T = $49.00, r = 3%, a T = 0.3, and T — 1.0 year were used. The curve for v = 40 
essentially overwrites the Black-Scholes curve for the scale used. The difference between the v — 40 and Black-Scholes 
prices were less then $0.10 for the range of So- 

Figure 9 gives the prices of calls and puts as a function of the t = price of the asset. The truncated Gosset 
formulae for Ct and Pt, Eqs. (9) and (10), were used to calculate the prices of the options. A comparison of Figs. 
8 and 9 shows that the truncated Gosset formula predicts lower prices for the options than does the capped Gosset 
formula. The truncated Gosset formula ignores events in the tails and thus costs less than the capped Gosset formula, 
which allows for events in the tails but places a limit on the value of the events in the tails. 
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FIG. 6: Price difference, Gosset — Black-Sholes, of a call 
option as a function of the confidence level for v = 3, 5, 
and 40 for a capped distribution. 
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FIG. 7: Price difference, Gosset — Black-Sholes, of a call 
option as a function of the confidence level for v = 3, 5, 
and 40 for a truncated distribution. 




FIG. 8: Prices of put and call options for v = 3 and v = 5 
using a capped Gosset formula with p — 0.999, and for 
the Black-Scholes formula. 




FIG. 9: Prices of put and call options for v = 3 and v = 5 
using a truncated Gosset formula with p — 0.999, and for 
the Black-Scholes formula. 



VII. CONCLUSION 



In this paper we have presented a method to price European options using a Student's ^-distribution. We refer to 
this method as a Gosset method and the formulae for pricing options as Gosset formulae, in honour of the person who 
published under the nom de plume (or pseudonym) Student. The Gosset formula essentially extends the Black-Scholes 
formula by removing the constraint that the volatility be constant in time. 

Student's ^-distributions have "fat tails" and fit the returns from assets better than a normal pdf. The difficulty 
with pricing options using the Student's £-distribution is that the fat tails cause one of the integrals, which is needed 
to price an option, to diverge. A viable and physically reasonable value for the price of the option can be obtained 
by either capping the value of the underlying asset or by truncating the distribution. The value for the cap or where 
the truncation takes place need not be unrealistic. The required integral can be evaluated for a confidence level of 
99.99%, which implies a price increase of > 4 000% for realistic parameters (volatility ctt = 0.4 and shape parameter 
v = 5) of the Student's ^-distribution. An increase of 4 000% in the value of an asset would be phenomenal and 
perhaps unphysical. 

There is little difference in philosophy between the approaches of pricing an option by capping the value of the 
asset or truncating the pdf. Both approaches reduce the cost of the option. However, there can be a large difference 
in the exposure to risk between the two approaches. 

Truncation of the pdf assumes that the value of the underlying asset will never exceed the value where the truncation 
occurs. There is no guarantee that this value will not be exceeded. 
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Capping the value of the asset ensures, by design and with consent, that the risk associated with the value of the 
option exceeding the value of the cap does not exist. In the cap approach, it is accepted that the value of the asset 
might exceed the value of the cap. However, the price of the option is set assuming that the payout will not exceed the 
value of the cap. Provided that the option writer and buyer agree to this approach, then there is no risk associated 
with a large increase in value of the asset. The buyer gives away the right to extreme returns in exchange for a lower 
price for the option. 

If the option writer and buyer do not agree to a cap on the payout, then the approach of capping the value of the 
asset still holds risk that the value of the asset will exceed the cap value. For the same critical value and parameters 
for the pdf, the cap approach offers an increased premium for the risk as compared to the approach of truncation of 
the pdf. This is because the cap approach allows for the possibility that the value of the asset will exceed the value 
of the cap. 

The Student's i-distribution has an additional parameter (the number of degrees of freedom, v) as compared to 
the normal pdf. This parameter is a statement of the belief about the range of values that the volatility might take 
at the time that the option expires. A large value of v is a statement that the volatility will change little in time. 
In the limit as v approaches infinity, the Student's i-distribution approaches the normal distribution and hence the 
Gosset formula for the price of an option becomes the same as the Black-Scholes formula. The Black-Scholes formula 
is derived on the assumption that the volatility is constant in time and known. A small value of v is a statement 
that the volatility might take one of a broad range of values. Fits to the DJIA for returns over 100 years give a 
value of v = 2.65. This value of v = 2.65 allows for the distribution of volatility that has been observed over the 100 
years of returns. In turbulent markets, it would be wise to use a small value for the number of degrees of freedom 
and a large confidence level. A confidence level is required in the pricing of options using a Student's i-distribution. 
In calm markets, it might be acceptable to use a large v, of order N — 1 where N is the number of independent 
samples that are used to calculate the volatility, and a smaller confidence level to match the expectations of the future 
and tolerance to risk. The flexibility offered by the Gosset formulae of choice of approach (cap or truncate) , of the 
number of degrees of freedom to match the fat tails of returns, and of the critical value allow the analyst and investor 
additional opportunity to use their insight and circumstances to price options. 

The approach used to derive the Gosset formulae can be adapted to distributions other than the Student's t- 
distribution. This may be of interest in situations and for returns on assets that do not appear to be well described 
by a Student's ^distribution. 
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